Draft version December 5, 2012 

Preprint typeset using L£-T^X style cmulatcapj v. 11/26/03 



(N 

o 

(N 
O 

Q 

6 



> 

(N 
OV 

o; 

(N 



X 



A NEW RADIO RECOMBINATION LINE MASER OBJECT TOWARD THE MONR2 HII REGION 
I. Jimenez-Serra 1 , A. Baez-Rubio 2 , V. M. Rivilla 2 , J. Martin-Pintado 2 , Q. Zhang 1 , M. Dierickx 1 and N. Patel 1 

Draft version December 5, 2012 

ABSTRACT 

We report the detection of a new radio recombination line (RRL) maser object toward the IRS2 source 
in the MonR2 ultracompact HII region. The continuum emission at 1.3 mm and 0.85 mm and the H30a 
and H26a lines were observed with the Submillimeter Array (SMA) at angular resolutions of ~0.5"-3" . 
The SMA observations show that the MonR2-IRS2 source is very compact and remains unresolved 
at spatial scales <400 AU. Its continuum power spectrum at millimeter wavelengths is almost flat 
(a=-0.16, with S v (x.v a ), indicating that this source is dominated by optically thin free-free emission. 
The H30a and H26a RRL emission is also compact and peaks toward the position of the MonR2-IRS2 
source. The measured RRL profiles are double-peaked with the H26a line showing a clear asymmetry 
in its spectrum. Since the derived line-to-continuum flux ratios (—80 and 180kms _1 for H30a and 
H26a, respectively) exceed the LTE predictions, the RRLs toward MonR2-IRS2 are affected by maser 
amplification. The amplification factors are however smaller than those found toward the emission 
line star MWC349A, indicating that MonR2-IRS2 is a weakly amplified maser. Radiative transfer 
modelling of the RRL emission toward this source shows that the RRL masers arise from a dense 
and collimated jet embedded in a cylindrical ionized wind, oriented nearly along the direction of the 
line-of-sight. High-angular resolution observations at sub-millimeter wavelengths are needed to unveil 
weakly amplified RRL masers in very young massive stars. 

Subject headings: stars: formation — masers — ISM: individual (Mon R2) — ISM: jets and outflows 



1. INTRODUCTION 

Calculations of the level populations of atomic hydro- 
gen in HII regions have shown that global population 
inversions can exist across the Rydberg levels of the hy- 
drogen atom 3 (Cillie 1936; Baker & Menzel 1938). Pop- 
ulation inversions lead to the formation of recombination 
lines under non-LTE conditions so that stimulated emis- 
sion can occur. The general understanding is that most 
HII regions in our Galaxy are optically thin at wave- 
lengths > 1 cm and stimulated emission is practically neg- 
ligible (Walmsley 1990). However, toward some ultra- 
compact (UC) HII regions, the presence of ionized stellar 
winds modifies the internal electron density structure of 
these sources (with n e >10 7 cm -3 ), leading to the forma- 
tion of optically thick cores where radio recombination 
line (RRL) masers can form (Martin-Pintado et al. 1989, 
1994). RRL maser amplification is therefore expected 
to be a common phenomenon in UC HII regions with 
evidence of stellar winds (Martin-Pintado 2002). 

The first RRL maser object was discovered by Martin- 
Pintado et al. (1989) toward the emission line star 
MWC349A. This source has an edge-on disk and a bipo- 
lar ionized flow (Cohen et al. 1985; Planesas et al. 1992; 
Weintroub et al. 2008; Martin-Pintado et al. 2011), whose 
continuum power law spectrum (a~0.6, with S v (xv a ; Al- 
tenhoff et al. 1981) is characteristic of a constant velocity 
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stellar wind (Olnon 1975). Martin-Pintado et al. (1989) 
showed that while the RRLs at 3 mm are faint with single 
Gaussian profiles, the RRLs at 1.3 mm are double-peaked 
with intensities factors of >50 brighter than the lines at 
3 mm. Since all RRLs with A<2 mm are largely amplified 
(Thum et al. 1998), Martin-Pintado et al. (1989) con- 
cluded that the RRLs toward MWC349A are masers. It 
has been proposed that r\ Carinae and Cepheus A HW2 
also show RRL maser emission (Cox et al. 1995; Jimenez- 
Serra et al. 2011). However, MWC349A still remains as 
the only RRL maser object firmly detected to date. 

In this Letter, we report the detection of a new RRL 
maser object toward the Monoceros R2 (MonR2) UC HII 
region with the SMA (d~830pc; Herbst & Racine 1976). 
MonR2 is a blister-type HII region (diameter of —27"; 
Massi et al. 1985; Wood & Churchwell 1989) that hosts a 
cluster of IR sources (e.g. Carpenter et al. 1997). Among 
these sources, IRS2 is a compact Young Stellar Object 
(YSO; Alvarez et al. 2004) with a luminosity of -5000 L© 
(Howard et al. 1994). This source (hereafter, MonR2- 
IRS2) is responsible for the spherical reflection nebula 
reported by Aspin & Walther (1990) in the MonR2 UC 
HII region. The detection of strong blue-shifted asymme- 
tries in the H26a RRL toward this source by the SMA, 
reveals that RRL maser amplification can form in dense 
ionized winds toward very young massive stars. 

2. OBSERVATIONS 

Observations of the H30a line (231.9 GHz) toward the 
MonR2-IRS2 source were carried out with the SMA 4 in 
the very extended (VEX) configuration in two tracks in 
single-receiver mode (4 GHz bandwidth per sideband). 

4 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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Table 1. Instrumental parameters of the SMA observations. 



Date 


Config. 


Line 


LO Freq. 


Synthesized Beam 


T225GHz 


±sys 


BP Cal. 


Flux Cal. 


Gain Cal. 








(GHz) 


("x", P.A.) 




(K) 








2010 Feb 13 


VEX 


H30a 


224.611 


0.53"x0.37", 48° 


0.05 


200-240 


3C273 


Titan 


0607-085/0730-116 


2010 Feb 20 


VEX 


H30a 


224.611 


0.53"x0.37", 48° 


0.03 


200-240 


3C273 


Titan 


0607-085/0730-116 


2011 Nov 15 


COM 


H30a 


226.227 


2.75" x 2.73", -77° 


0.07 


200-240 


BLLAC 


Ganymede 


0607-085/0530+135 


2011 Nov 15 


COM 


H26a 


348.324 


2.29" x 1.50", -39° 


0.07 


400-500 


BLLAC 


Ganymede 


0607-085/0530+135 



In addition, the H30a and the H26a (353.6 GHz) RRLs 
were simultaneously observed in a third track in com- 
pact (COM) configuration in dual-receiver mode, which 
provided 2 GHz of total bandwidth per sideband and re- 
ceiver. The instrumental parameters of the SMA ob- 
servations are reported in Table 1. The phase center 
of the observations was set at a(J2000)=06 /l 07 m 45.83 s , 
<5(J2000)=-06°22'53.50". We used a uniform spectral 
resolution of 0.8 MHz, which provided a velocity resolu- 
tion of -1.1 km s -1 at 231.9 GHz, and of -0.7 km s" 1 at 
353.6 GHz. Data calibration was carried out within the 
IDL MIR software package, and continuum subtraction, 
imaging and deconvolution was done within MIRIAD. 
The uncertainty in the flux calibration was within 20%. 

3. RESULTS 

In Figure 1 (left panels) , we present the continuum im- 
ages at 1.3 mm measured toward MonR2-IRS2 in COM 
and VEX configurations (grey scale; upper and middle 
panels), and at 0.85mm in COM (lower panel). The de- 
rived parameters of this emission, obtained by perform- 
ing 2D Gaussian fits, are given in Table 2. The position 
of the continuum peak is very close to that reported for 
the IRS2 source in the near-IR and X-rays (Carpenter 
et al. 1997; Nakajima et al. 2003). The derived angular 
sizes for the MonR2-IRS2 source indicate that this ob- 
ject is very compact and unresolved in the SMA images 
(Table 2). Since the continuum peak intensity at 1.3 mm 
is -0.15 Jy beam^ 1 in both VEX and COM, MonR2- 
IRS2 does not show any structure at sub-arcsecond scales 
(0.4" -0.5" or 320-400 AU). This is consistent with the re- 
sults by Alvarez et al. (2004) from near-IR speckle imag- 
ing. By smoothing the 1.3 mm and 0.85 mm images in 
COM to the same angular resolution (—2.8"), we derive 
a decreasing spectral index of a=— 0.16 (with S u (xv a ) for 
MonR2-IRS2. This spectral index is consistent with op- 
tically thin free-free continuum emission. We note that 
our estimate of the spectral index is not affected by miss- 
ing flux because the MonR2-IRS2 source is very compact. 
Indeed, its measured flux at 1.3 mm is the same in the 
VEX and COM beams (see Table 2), confirming the lack 
of large-scale structures in this source. 

Superimposed on the continuum images, Figure 1 re- 
ports the integrated intensity emission of the H30a 
and H26a lines measured toward MonR2-IRS2 from 
-35.5 kms -1 to 47.0kms~ 1 (thick contours). The RRL 
maps show a compact structure centered at the contin- 
uum peak, which is almost unresolved even in the VEX 
images at angular scales of —0.5". The spectra of the 
H30a and H26a RRLs extracted from the position of the 
continuum peak are shown in the right panels of Figure 1. 
These spectra have been smoothed to a velocity resolu- 
tion of 2.2kms _1 for H30a, and of 2.8 kms" 1 for H26a. 
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Fig. 1. — Left panels: Integrated intensity images of the H30q 
line from -35.5 to 47.0kms _1 (thick contours) observed with the 
SMA in COM and VEX (upper and middle panels), and of the 
H26a RRL measured toward MonR2-IRS2 in COM (lower panel). 
These images are superimposed on the continuum emission ob- 
served at 1.3 mm in COM and VEX (grey scale; upper and middle 
panels) and at 0.85 mm in COM (lower panel). Offsets are given 
in arcseconds with respect to the phase center of the observations. 
Note that the spatial scales shown for the VEX images (middle 
panel) are ~5 times smaller than those of the COM data (upper 
and lower panels). The first contour and step level for the 1.3 mm 
and H30a maps in COM (upper panel) are, respectively, 13.5 (3c) 
and 27 mjy beam -1 , and 2.4 (3<r) and 2.4 Jy beam -1 km s" 1 . For 
the 1.3 mm and H30ct maps in VEX (middle panel), the first con- 
tour and step level are 4.5 (3cr) and 18mjybeam , and 0.75 (3c) 
and 2.25 Jy beam -1 kms -1 , respectively. For the 0.85 mm and 
H26a images in COM (lower panel), these levels are 30 (3<r) and 
30mjybeam _1 , and 4.5 (3<x) and 4.5 Jy beam -1 kms~ 1 , respec- 
tively. Right panels: Spectra of the H30a and H26a RRLs mea- 
sured with the SMA in COM and VEX toward the MonR2-IRS2 
continuum peak. Dashed curves show the individual Gaussian fits 
of the different velocity components detected toward MonR2-IRS2. 
Solid lines show the total Gaussian fit of the H30a and H26a RRLs. 
The vertical dotted line indicates the radial velocity of the source 
(?, iSH =10.5kms- 1 ; Torrefies et al. 1983). 
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Table 2. Derived parameters of the continuum and RRL emission to- 
ward MonR2-IRS2. 







Continuum Emission 




Config. 


A 


Position Cont. 


Peak 


Peak Flux a 


Angular Size 




(mm) 


a(J2000) 


5(J2000) 


(Jy beam -1 ) 


("x", PA.) 


OW1VI 


1 Q 




-06°22'53.50" 


0.154±0.005 


o.lo X z. /y , lUo 


VEX 


1.3 


06 h 07 m 45 .806 s 


-06°22'53.53" 


0.157±0.002 


0.55" x 0.38", 45° 


COM 


0.85 


06 h 07 m 45.807 s 


-06°22'53.45" 


0.14±0.01 


2.44" x 1.60", 130° 








RRL Emission 




Config. 


Line 


T L Av h 


VLSR 


Ad 


T L C 






(Jy beam - 1 kms -1 ) 


(kms - ) 


(kms -1 ) 


(Jy beam x ) 


COM 


H30a 


6.93 (0.22) 


-13.4 (0.6) 


23.4 (1.7) 


0.28 (0.03) 






8.36 (0.24) 


10.5 (0.0) 


30.0 (0.0) 


0.26 (0.03) 






4.69 (0.19) 


33.3 (0.6) 


18.9 (2.0) 


0.23 (0.03) 


VEX 


H30a 


4.2 (0.4) 


-10.89 (0.12) 


15.4 (0.3) 


0.26 (0.07) 






4.7 (0.6) 


10.5 (0.0) 


30.0 (0.0) 


0.15 (0.07) 






3.3 (0.4) 


31.46 (0.15) 


13.7 (0.4) 


0.23 (0.07) 


COM 


H26a 


14.5 (0.5) 


-14.0 (3.0) 


15.5 (3.0) 


0.88 (0.05) 






6.7 (0.7) 


10.5 (0.0) 


32.2 (3.0) 


0.20 (0.05) 






5.0 (0.4) 


28.9 (3.0) 


10.4 (3.0) 


0.45 (0.05) 



a The error in the peak flux corresponds to the Icr noise level in the continuum images. 

b The error in the integrated intensity flux of the RRLs is calculated as a ^^ rea = lcrXV 'Sv X Av, with Sv the velocity resolution in the RRL 
spectra (see Section 3), and Av the linewidth derived from the Gaussian fit of the RRL emission. 
c The error in the peak intensity of the H30a and H26ct lines is given by the lcr noise level in the RRL images. 



The H30ct and H26a lines have bright double-peaked line 
profiles whose components are red- and blue-shifted by 
^20-25 kms -1 with respect to the radial velocity of the 
source (wls,r^10.5 kms -1 ; Torrelles et al. 1983). In ad- 
dition, some faint emission at WLSfl~10-5kms -1 is also 
detected, likely associated with the extended MonR2 UC 
HII region. This would explain why the double-peaked 
line profile of H30o; is not as clearly seen in COM as in 
VEX, since the contribution from the extended emission 
would be larger within the COM beam. 

From Figure 1 , we also find that while the double- 
peaked H30a RRLs have similar peak intensities of ^0.2- 
0.3 Jy beam -1 , the H26a line profile shows a clear asym- 
metry with the blue-shifted gas being factors of ^2 
brighter than the red-shifted emission (^0.9 Jy beam -1 
vs. ^0.5 Jy beam -1 , respectively). This behavior can 
only be explained if stimulated emission plays a key role 
in the formation of the H26a RRL (see Martin-Pintado 
et al. 1993; Jimenez-Serra et al. 2011; Baez-Rubio et al. 
2012). As shown in Section 5, the H26a RRL blue-shifted 
asymmetry is a consequence of the presence of more back- 
ground emission to be amplified by the blue-shifted fore- 
ground ionized gas than for the red-shifted material, as- 
suming that the RRL masers are formed in an expanding 
ionized jet. 

The observed parameters of the three velocity com- 
ponents of the H30a and H26a RRLs are reported in 
Table 2. These parameters were derived by fixing the 
peak velocity and linewidth of the ionized component at 
ambient velocities to, respectively, 10.5 kms -1 (the ra- 
dial velocity of the MonR2-IRS2 source) and to 30 km s -1 
(the expected thermal linewidth for RRLs in HII regions 
with electron temperatures of X" e =l-2xl0 4 K; see e.g. 
Keto et al. 2008). The derived peak velocities of the 
H30a and H26a components are 29-33 kms -1 for the 
red-shifted gas, and —11 to —14 kms -1 for the blue- 
shifted emission. The measured linewidths are ^10- 
15 kms -1 (Table 2), except for the H30a COM data 



whose linewidths are broader due to the contamination 
from the extended MonR2 UC HII region. 

By integrating the emission of the H30a RRL from -18 
to -4.8 kms -1 and from 27.1 to 38.1 kms -1 , the VEX im- 
ages of the blue-shifted and red-shifted H30a gas reveal 
a spatial shift of 0.045" (or 360 AU) in the southeast- 
northwest direction (blue-shifted emission toward the 
southeast, red-shifted gas toward the northwest). Al- 
though this shift is small compared to the angular reso- 
lution of the VEX data (0.5" x0.4"; Table 1), higher po- 
sitional accuracy of 2xS/n ( wl th S/N the signal-to- noise 
ratio of the H30a emission) can be achieved thanks to the 
bright integrated intensities of the blue- and red-shifted 
H30a RRL components. By fitting 2D Gaussians to the 
integrated intensity images of these velocity components, 
we obtain an accuracy in the H30a Gaussian centroid po- 
sition of ^0.008". This implies that the detected H30a 
spatial shift of 0.045" ±0.008" is at the 5.5cr confidence 
level. This shift could be associated with a collimated 
ionized jet propagating at a velocity of 20 kms -1 (see 
Section 5). 

4. RRL MASERS IN MONR2-IRS2 

The double-peaked RRL profiles with clear asymme- 
tries in the low-n RRL transitions toward MonR2-IRS2 
(i.e. in the H26a line; Section 3), resemble the be- 
haviour of the RRLs detected toward MWC349A, and 
suggest that the RRL emission toward MonR2-IRS2 
at A<1.3mm is affected by maser amplification. The 
integrated line-to-continuum flux ratios, ILTRs 5 , de- 
rived toward this source support this idea (Table 3). 
As in MWC349A (Martin-Pintado et al. 1989; Thum 
et al. 1994), the ILTRs derived toward MonR2-IRS2 
(80 kms -1 for H30a and 180 kms -1 for H26a) exceed the 

5 ILTR is defined as AvT^/Tq, and depends on the RRL fre- 
quency as v 1 1 for LTE and optically thin emission (see e.g. Martin- 
Pintado et al. 1989). 
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Table 3. Integrated Line-to-Continuum flux ratios (ILTRs). 



RRL 


LTE a 


MonR2-IRS2 b 


MWC349A 




(kms -1 ) 


(km s —1 ) 


(kms- 1 ) 


H3GV* 


~65 


80 


298 c 


H26a 


~103 


180 


756 d 



a Calculatcd for optically thin continuum emission, T* = 10 4 K and 
AT(Hc+)/7V(H+)=0.08. 

b Estimated from the ratio between the total area of the H30a and 
H26a lines and the continuum peak flux at 1.3 mm and 0.85 mm 
(Table 2). 

c From Martin-Pintado et al. (1989). 

d Dcrivcd from the average H26a line flux measured by Thum et al. 
(1994), and from the expected frcc-frcc continuum flux at 0.85 mm 
assuming a spectral index of a~0.6 (Altcnhoff et al. 1981). 
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Fig. 2. — Comparison of the H30ct and H26a RRL spectra mea- 
sured with the SMA toward MonR2-IRS2 (in VEX and COM, re- 
spectively), with those predicted by the MORELI code for two 
geometries of the MonR2-IRS2 source: a) a cylindrical, ionized 
outflow oriented along the line-of-sight (left panels); and b) a cylin- 
drical, ionized outflow with an elongated and dense inhomogeneity 
resembling a collimated jet (right panels). Vertical dotted lines in- 
dicate the assumed radial velocity of the source t)_LSi?=10kms _1 . 



values predicted under LTE conditions (65kms -1 and 
103 km s -1 , respectively). However, unlike MWC349A, 
the RRL maser amplification factors toward MonR2- 
IRS2 with respect to LTE are significantly smaller than 
those found in MWC349A (-1.2-1.7 for MonR2-IRS2 vs. 
>5-7 for MWC349A; Table 3 and Strelnitski et al. 1996; 
Thum et al. 1998). This indicates that MonR2-IRS2 is a 
weakly amplified RRL maser. This is due to the fact that 
MWC349A has a larger optically thick core than MonR2- 
IRS2, allowing a much larger maser amplification of the 
RRLs (see Ponomarev et al. 1994, and Section 5). 

5. MODELLING OF THE RRL EMISSION TOWARD 
MONR2-IRS2 

By comparing the H30a and H26a RRLs measured to- 
ward MonR2-IRS2 with radiative transfer modelling of 
the RRL emission, we can constrain the physical struc- 
ture and kinematics of the ionized gas toward this YSO. 
To do this, we have used the three-dimensional RRL 
radiative transfer code MORELI developed by Marti'n- 
Pintado et al. (2011) and Baez-Rubio et al. (2012). In our 
calculations, we have used the departure coefficients, b n 
and /3„, derived by Walmsley (1990). For the geometry 
of the source, we have assumed that MonR2-IRS2 has a 
central mass of 13 M©, which is consistent with a Bl star 
on the ZAMS (luminosity of -5000 L Q ; Panagia 1973) 



and with the Lya photon flux (~ 6 x 10 45 photons s 1 ) 
derived from the 1.3 mm and 0.85 mm continuum data 
(Table 2). In addition, we consider that MonR2-IRS2 
powers a cylindrical, isothermal (T e =10 4 K) ionized out- 
flow, which expands at a velocity of 20 km s _1 . The cylin- 
drical outflow is defined so that its outer radius is 100 AU 
and its electron density distribution n e decreases as r~ 2 , 
with r the distance to the star and with n e =10 8 cm~ 3 
at a distance r— 7AU. The total length of the cylinder 
is L=520 AU and its axis is oriented along the direction 
of the line-of-sight. The assumed radial velocity for the 
MonR2-IRS2 source is VLSit^lOkms -1 . 

To reproduce the optically thin free-free continuum 
emission of MonR2-IRS2 at 1.3 mm and 0.85 mm, we 
consider a cavity around the star with radius r=13AU 
(see Baez-Rubio et al. 2012). The predicted continuum 
fluxes at 1.3 mm and 0.85 mm are 0.156 Jy beam" 1 and 
0.148 Jy beam -1 , respectively, which differ from the ob- 
served values by less than 3%. 

In Figure 2 (left panels) , we compare the RRL profiles 
of the H30a and H26a lines measured with the SMA 
toward the MonR2-IRS2 source (in VEX and COM, re- 
spectively), with the synthetic line profiles predicted by 
the MORELI code for the cylindrical, ionized outflow as- 
sumed for MonR2-IRS2 (Case a). Although our model 
reproduces the double-peaked profile of the H30a RRL 
relatively well, it clearly fails to predict the asymmetry 
detected in the H26a line (Figure 2). Indeed, the pre- 
dicted line profiles for both the H30a and H26a RRLs are 
symmetric because the continuum emission of MonR2- 
IRS2 is optically thin throughout the ionized outflow, 
which prevents the strong amplification and the asym- 
metries of the RRLs. 

In order to explain the strongly asymmetric H26a line 
profile, larger electron densities are required for the in- 
ner regions of the ionized wind so that the continuum 
emission becomes locally optically thick. To simulate 
this, we have assumed two elongated inhomogeneities 
with n e =10 8 cm -3 located at the axis of the blue-shifted 
and red-shifted lobes of the cylindrical ionized outflow, 
resembling a collimated ionized jet. The size of the inho- 
mogeneities is 2 AUx2 AUx60 AU, and they are placed 
on the axis of the cylindrical outflow at a distance of 
43 AU from the central star. 

In Figure 2 (right panels), we report the H30a and 
H26a line profiles predicted by our model with a cylin- 
drical ionized outflow powered by a dense, collimated jet 
(Case b). While the predicted H30a emission shows a 
similar line profile to that obtained in Case a) for only 
a cylindrical ionized outflow, the model with the ionized 
outflow+jet perfectly matches the asymmetry observed 
in the H26a line profile. This is due to the fact that the 
inner regions of the jet are optically thick, allowing the 
maser amplification of the H26a line. Since the optically 
thick region in the plane of the sky of the ionized jet is 
very small (only 2 AUx2 AU), its contribution to the to- 
tal optical depth of the MonR2-IRS2 continuum source is 
practically negligible. We note that only the H26a RRL 
gets amplified because the departure coefficient (3 n of the 
H26a transition at an electron density of n e =10 8 cm~ 3 
is significantly larger than that of the H30a line (see 
Figure 8 in Strelnitski et al. 1996). The electron den- 
sities n e assumed for the elongated inhomogeneities do 
not likely exceed 10 8 cm -3 , because the /3„ coefficients 
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for both the H30a and H26cn RRLs get close to zero, 
hindering the RRL maser amplification (see Figure 5 in 
Strelnitski et al. 1996). 

In our model for MonR2-IRS2 with an ionized out- 
flow +jet (Case b), the elongation of the dense inho- 
mogeneities is required to sufficiently amplify the blue- 
shifted peak of the H26a emission. In contrast with the 
blue-shifted peak, the RRL emission from the red-shifted 
lobe of the ionized jet does not get amplified because the 
optically thick continuum acts as a screen for the red- 
shifted H26a emission. 

Finally, we stress that other geometries for the MonR2- 
IRS2 source have been explored in our study. In the case 
of an edge-on Kcplerian rotating disk with a bi-conical 
ionized flow (similar to that assumed for MWC349A; 
see Martin-Pintado et al. 2011), the model also predicts 
double-peaked line profiles for large opening angles of 
the ionized wind (# a >60°; sec Baez-Rubio et al. 2012, 
for the definition of this angle). However, the model not 
only fails to reproduce the blue-shifted asymmetry of the 
H26a line profile, but predicts a RRL emission excess at 
velocities close to the ambient cloud velocity that is not 
observed. The presence of high-density inhomogencitics 
could partially alleviate the discrepancies between the 
model predictions and the observations, but the best fit 
is obtained with the cylindrical ionized outflow powered 
by the elongated jet. 

In summary, we report the detection of a new RRL 
maser object toward the IRS2 source in the MonR2 UC 
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at sub-millimeter wavelengths have the potential to un- 
veil a population of weakly amplified RRL maser objects 
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